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ABSTRACT 

Direct and indirect pieces of observational evidence point to a strong connection be- 
tween high-redshift quasars and their host galaxies. In the framework of a model 
where the shining of the quasar is the episode that stops the formation of the galactic 
spheroid inside a virialized halo, it has been proven possible to explain the submil- 
limetrc source counts together with their related statistics and the local luminosity 
function of spheroidal galaxies. The time delay between the virialization and the quasar 
manifestation required to fit the counts is short and incresing with decresing the host 
galaxy mass. In this paper we compute the detailed chemical evolution of gas and stars 
inside virialized haloes in the framework of the same model, taking into account the 
combined effects of cooling and stellar feedback. Under the assumption of negligible 
angular momentum, we are able to reproduce the main observed chemical properties 
of local ellipticals. In particular, by using the same duration of the bursts which are 
required in order to fit the submillimetre source counts, we recover the observed in- 
crease of the Mg/Fe ratio with galactic mass. Since for the most massive objects the 
assumed duration of the burst is T^urst ^ 0.6 Gyr, we end up with a picture for el- 
liptical galaxy formation in which massive spheroids complete their assembly at early 
times, thus resembling a monolithic collapse, whereas smaller galaxies are allowed for 
a more prolonged star formation, thus allowing for a more complicated evolutionary 
history. In the framework of the adopted scenario, only quasar activity can provide 
energies large enough to stop the star formation very soon after virialization in the 
most massive galactic haloes. The chemical abundance of the gas that we estimate 
at the end of the burst matches well the metallicity inferred from the quasar spectra. 
Therefore, the assumption that quasar activity interrupts the main episode of star 
formation in elliptical galaxies turns out to be quite reasonable. In this scenario, we 
also point out that non-dusty extremely red objects are the best targets for searching 
for high-redshift Type la supernovae. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution ~ galaxies: for- 
mation - galaxies: ISM - galaxies: stellar content - nuclear reactions, nucleosynthesis, 
abundances 



1 INTRODUCTION 

A traditional view for the formation history of an ellipti- 
cal galaxy stems from the Milky Way collapse model of 
Eggen, Lynden-Bell & Sandage (1962): monolithic collapse 
and rapid star formation lead to a subsequent track known 
as 'passive evolution' (i.e., without further star formation). 
This behaviour can be obtained if a 'galactic wind' develops 
as a result of the energy injection into the gas by supernovae 
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(SNe): the gas still present in the galaxy is swept away and 
the subsequent evolution is determined only by the amount 
of matter and energy which is restored into the interstellar 
medium (ISM) by the dying stars (Mathews & Baker 1971; 
Larson 1974a; Ikeuchi 1977). The characteristic time-scales 
for the occurrence of the galactic winds depend on the pre- 
vious star formation history, on the assumptions on the dark 
matter amount and distribution, and on the SN energetics. 
This scenario is known as early monolithic collapse (e.g., 
Larson 1974b; Larson & Tinsley 1974; Arimoto & Yoshii 
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1987; Matteucci & Tornambe 1987; Matteucci 1992, 1994; 
Bressan, Chiosi & Fagotto 1994). 

On the other hand, hierarchical clustering in the Uni- 
verse is a fundamental process. The merging events, taking 
place over a major fraction of the cosmological time, have 
been associated with bursts of star formation in building-up 
galaxies in semi-analytical models (e.g., Kauffmann, Guider- 
doni & White 1994; Baugh, Cole & Frenk 1996; Kauffmann 
1996; Kauffmann & Chariot 1998). 

A long-lasting star formation has been thoroughly ruled 
out, at least in giant ellipticals, on the basis of the chemical 
properties of the stellar population and the gas leftover by 
the process of star formation, which constrain the chemical 
enrichment history of these spheroids, and hence the mecha- 
nisms responsible for their formation (Bower, Lucey & Ellis 
1992; Ellis et al. 1997; Bernardi et al. 1998; Thomas 1999; 
Thomas & Kauffmann 1999). Studies on luminosity func- 
tions at substantial redshift show that the majority of lumi- 
nous E and SO galaxies were already in place at z ~ 1 and 
that most of their stars were relatively old (Im et al. 2001; 
Cohen 2001). Searches for a population of faint intrinsically 
red objects, representing the expected z > 1 precursors of 
passively evolving ellipticals which formed their stars at high 
redshift, have been conducted: Daddi et al. (2000) find that 
colours and surface density of extremely red objects (EROs) 
are consistent with the hypothesis that they are the precur- 
sors of the present day luminous elliptical galaxies. Also the 
clustering properties support this conclusion (Daddi et al. 
2000; McCarthy et al. 2000; Magliocchetti et al. 2001). Sev- 
eral authors (e.g., Franceschini et al. 1991, 1994; Lilly et al. 
1999; Dunlop 2001; Granato et al. 2001) suggested that the 
burst, during which most of the stars in spheroidal galaxies 
were formed at 2 > 1, is elusive in the optical bands, but 
shows up in the far-lR domain. Evidence of this phase has 
been found by Small, Ivison & Blain (1997) as a result of a 
deep survey with SCUBA and has been confirmed by sub- 
sequent surveys and foUows-up (Hughes et al. 1998; Barger 
et al. 1998; Blain et al. 1999; Small et al. 2000; Almaini et 
al. 2001; Borys et al. 2001). 

The above-sketched picture is complemented by addi- 
tional pieces of observational evidence pointing to a strong 
connection between quasi-stellar objects (QSOs) and galax- 
ies. In particular, the recently discovered correlation be- 
tween the central massive dark object (MDO) (routinely 
interpreted as a dormant black hole) and the hot stellar 
component of nearby galaxies (e.g., Kormendy & Richstone 
1995; Magorrian et al. 1998; van der Marel 1999; Ferrarese 
& Merritt 2000; Gebhardt et al. 2000; McLure & Dunlop 
2001; Merritt & Ferrarese 2001a, b) suggests a direct connec- 
tion between QSO activity and galaxy formation. Since ~ 
97 per cent of early-type galaxies do have MDOs (Magorrian 
et al. 1998), it seems reasonable to construct models taking 
into account jointly the cosmological formation of QSOs and 
spheroids (Silk & Rees 1998; Friaga & Terlevich 1998; Kauff- 
mann & Haehnelt 2000; Monaco, Salucci & Danese 2000; 
Granato et al. 2001, among others). While QSO activity and 
evolution at low and intermediate redshifts are thought to 
be driven by interactions between galaxies (Barnes & Hern- 
quist 1991), the QSO evolution at z » 1 is expected to 
be more closely related to the formation of massive galaxies 
in overdense regions (Haehnelt & Rees 1993). In particu- 
lar, Granato et al. (2001) showed that most of the above- 



mentioned observations can be explained in a hierarchical 
model for spheroidal galaxy formation, if stellar and QSO 
feedback are properly included. In particular, they pointed 
out that the stellar feedback is more effective in reducing 
the star formation rate of smaller haloes, thus producing a 
star formation rate which is increasing with the halo mass 
and the final mass in stars. 

The present study aims to investigate further the evolu- 
tionary histories for elliptical galaxies within the framework 
suggested by Granato et al. (2001), focusing on the expected 
chemo-photometric properties of the resulting stellar popu- 
lations and on the chemistry of the gas. The predictions of 
the model are then tested against the body of the available 
data. We assume that after a vigorous star formation build- 
ing up the bulk of the stellar population and enriching the 
gaseous medium to roughly solar or higher metallicity, the 
QSO shines at the centre, ionizing the surrounding medium 
and inhibiting further star formation. A galactic wind is pos- 
sibly triggered at this point, due to the combined feedback of 
QSO and stars. The QSOs shine in an inverted hierarchical 
order (Monaco et al. 2000), which means that the time in- 
terval between the onset of the star formation and the peak 
of the QSO activity is shorter for more massive spheroids. 

The paper is organized as follows: in §2 we report on the 
current status of observations, as far as both the integrated 
properties of the composite stellar population and the inter- 
stellar/surrounding gaseous component are concerned. In §3 
we introduce our model and in §4 we present and compare 
the model results to the available observations. A summary 
is given in §5, together with a comparison with previous 
studies. In the following we will assume Q,a ~ 0.7, SIm ~ 0.3 
and h = 0.7, unless otherwise stated. 



2 OBSERVATIONAL CONSTRAINTS 
2.1 The stellar populations 

Most of the information on elliptical galaxies comes from 
their integrated properties: abundances are derived either 
through colours or integrated spectra. The metallicity in- 
dices Mg2 and (Fe), as originally defined in Faber, Burstein 
& Dressier (1977) and Faber et al. (1985), are the most 
widely used metallicity indicators and have been measured 
for a number of objects. The integrated properties of ellipti- 
cal galaxies can be analysed by means of population synthe- 
sis techniques (e.g., Buzzoni, Gariboldi & Mantegazza 1992; 
Bruzual & Chariot 1993; Bressan, Chiosi & Fagotto 1994; 
Bressan, Chiosi & Tantalo 1996; Tantalo, Chiosi & Bressan 
1998) in order to get an estimate of the real abundances. Pos- 
itive values of the elemental ratio [Mg/Fe] have been found in 
the nuclei of ellipticals, increasing with increasing the galac- 
tic mass (O'Connell 1976; Peterson 1976; Faber, Worthey 
& Gonzalez 1992; Worthey, Faber & Gonzalez 1992; Weiss, 
Peletier & Matteucci 1995; Tantalo et al. 1998; Worthey 
1998; J0rgensen 1999; Kuntschner 2000; Kuntschner et al. 
2001; Terlevich & Forbes 2001). Although large uncertain- 
ties are involved with these methods, the overabundance 
of the [Mg/Fe] ratio in ellipticals with respect to the so- 
lar value is accepted and generally interpreted as due to a 
short, intense star formation, perhaps coupled to an initial 
mass function (IMF) biased towards massive stars. Other 
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observational clues consistent with this picture are reminded 
below: 

Bender, Burstein & Faber (1992) studied the distribu- 
tion of elliptical galaxies in the Virgo cluster and in the 
Coma cluster in a 3-space of the basic global parameters 
central velocity dispersion a, effective radius r^, and mean 
effective surface brightness le. They concluded that clus- 
ter ellipticals lie on a fundamental plane (FP, Djorgovski & 
Davis 1987), i.e., they do not fill uniformly the space of the 
global parameters but rather select a bidimensional region 
of it (see also Dressier et al. 1987). The tightness of the FP 
for elliptical galaxies in local clusters was then confirmed 
by other works (e.g., Renzini & Ciotti 1993). ft points to a 
high uniformity and synchronization in the galaxy formation 
process in rich clusters. 

Bower et al. (1992) pointed out the existence of a very 
tight colour -cr relation for elliptical galaxies in the Virgo 
and Coma clusters. By combining these cluster data at low 
redshift with /fST-selected samples at intermediate redshift, 
it was demonstrated that at least cluster ellipticals are made 
of very old stars, with the bulk of them having formed at 
redshift z > 2 (Ellis et al. 1997). A tight colour - magnitude 
relation (CMR) was also found, in the sense that colours 
become redder with increasing luminosity. The tightness of 
the CMR for elliptical galaxies in clusters has been con- 
firmed up to z ~ 1 (Aragon-Salamanca et al. 1993; Stan- 
ford, Eisenhardt & Dickinson 1998). The optical-lR colours 
of the early-type cluster galaxies observed by Stanford et al. 
(1998) become bluer with increasing z (out to z — 0.9) in a 
manner consistent with the passive evolution of an old stel- 
lar population formed at an early cosmic epoch. Moreover, 
the modest shift with increasing z found in the zero-point of 
the FP, Mg2 - cr, and CMRs of cluster ellipticals (e.g., Dick- 
inson 1995; Ellis et al. 1997; Ziegler & Bender 1997; Bender 
et al. 1998; van Dokkum et al. 1998; Kodama et al. 1998; 
Stanford et al. 1998) has leaded to the conclusion that most 
stars in cluster ellipticals formed at z > 3. 

However, until a couple of years ago much of this evi- 
dence was restricted to cluster elliptical galaxies. Using new 
observations for a sample of 931 early-type galaxies assigned 
to three different environments (clusters, groups, and field), 
Bernardi et al. (1998) demonstrated that cluster, group, and 
field early- type galaxies follow almost identical Mg2 - a re- 
lations. The largest Mg2 zero-point difference they found is 
only 0.007 ± 0.002 mag, thus implying an age difference 
of only ~ 1 Gyr between the stellar populations of cluster 
and field elliptical galaxies, these latter being younger on 
the average. Recently, it has been claimed that significant 
age variations might be hidden in the Mg2-(T relation, so 
that it would be desirable to consider the two Balmer line 
indices H/3 and Hn/Fe as more robust age indicators (Con- 
cannon, Rose & Caldwell 2000). However, the scatter of H/3 
observed in cluster and luminous field elliptical galaxies can 
also be explained by composite populations that contain a 
small fraction of old metal-poor stars, without requiring a 
young stellar component (Maraston & Thomas 2000). On 
the other hand, there is evidence in both cluster and field 
ellipticals that a modest fraction of the present day stellar 
population is indeed formed in secondary bursts at relatively 
low redshift (e.g., Trager et al. 2000a,b; Treu et al. 2001). 

Another piece of observational evidence supporting a 
scenario where an early, intense star formation is followed 



by a longer, quiescent phase (passive evolution) is given by 
the present day Type la and H SN rates in early-type galax- 
ies: Rates N la = 0.18 ± 0.06 SNu|1] Rates nii < 0.02 SNu 
(for Ho — 75 km Mpc~^; Cappellaro, Evans & Turatto 
1999). Since Type H SNe come from short-living progeni- 
tors, whereas Type la SNe come from long-living ones, the 
observed Type la and H SN rates imply that the star forma- 
tion in early-type galaxies must be inactive at the present 
time. 

2.2 The gaseous medium 

2.2.1 The quasar environment 

In a recent paper, Hamann & Ferland (1999) have reviewed 
the topic of QSO emission and intrinsic absorption lines. 
They have shown that these lines can give us some hints 
on the chemical composition of the QSO environment. The 
metallicity of broad emission line regions (BELRs) can be 
studied by using emission ratios. Hamann et al. (2001) find 
evidence for roughly solar or higher metallicities out to z > 
4 and enhanced N/C ratios with respect to solar in more 
luminous objects (confirmed by Fan et al.'s 2001 spectral 
analysis for a few z ~ 6 quasars). This implies that most 
of the enrichment of the gaseous medium must occur before 
the QSOs become observable (i.e., on time-scales < 1 Gyr, 
at least at the highest redshifts). This gives further support 
to the idea that QSOs are hosted at the centres of elliptical 
galaxies that experienced vigorous star formation at some 
early stage of their evolution. 

2.2.2 The hot gas 

The standard SN rates predict the iron abundance of the 
ISM in early-type galaxies to be as high as several times 
the solar value (e.g., Loewenstein & Mathews 1991; Ciotti 
et al. 1991; Renzini et al. 1993). However, the first measure- 
ments of the ISM with ASCA showed that the metallicity 
was less than half a solar (Awaki et al. 1994; Loewenstein 
et al. 1994; Mushotzky et al. 1994; Matsushita et al. 1994; 
Arimoto et al. 1997; Matsumoto et al. 1997). Giant early- 
type galaxies are estimated to have a roughly solar stellar 
iron abundance (e.g., Arimoto et al. 1997). Thus, the X-ray 
measured abundances implied that the ISM metal abun- 
dance was even lower than the stellar metallicity. Lately, 
from ASCA and 7? 05^4 T spectra of four among the brightest 
elliptical galaxies in X-rays, Fe abundances of 1-2 times 
solar and (except for one object) relative abundances fixed 
at their solar values have been found (Buote 1999)^ This 
result has been confirmed by the analysis of ASCA data for 
27 giant early-type galaxies by Matsushita, Ohashi & Mak- 
ishima (2000): they find nearly solar (within a factor of 2) 
abundances of Fe and Q-elements in X-ray luminous galax- 
ies. Much more room for uncertainty is left for X-ray fainter 
galaxies, but it seems that the contribution from Type la 
SNe to the ISM abundance is lower in X-ray faint than in 

t 1 SNu = 1 SN per century per 10^° Lbq- 

t The Fe abundances obtained by Buote should be increased by a 
factor of 1.44 to reflect the meteoritic solar abundances instead of 
the photospheric values he used (see Ishimaru &: Arimoto 1997). 
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X-ray luminous systems. These results strongly suggest that 
a large fraction of the SNIa products have escaped into the 
intergalactic space. 

An additional clue to the possible relation between 
QSOs and galaxy formation comes from X-ray observations 
of the intergalactic medium (IGM) in galaxy clusters. As 
soon as deviations from the expected Lx oc relation 
emerged from the cluster data, several authors suggested 
a possible pre-heating of the intergalactic gas (e.g. Kaiser 
1991; Evrard & Henry 1991). Evidence of an entropy ex- 
cess with respect to the values expected from purely gravi- 
tational heating has been recently found by Ponman, Can- 
non & Navarro (1999) and Lloyd-Davies, Ponman & Can- 
non (2000), although the interpretation of the data is still 
debated (see, e.g., Roussel, Sadat & Blanchard 2000). The 
extra-energy injected is in the range ~ 0.5-3 keV per par- 
ticle (e.g., Cavaliere, Menci & Tozzi 1997; Wu, Fabian & 
Nulsen 2000; Borgani et al. 2001). On one hand, the most 
natural heating mechanism, heating from SNe, seems to fall 
short of the required energy budget (Valageas & Silk 1999; 
Kravtsov & Yepes 2000 and refs. therein; but see also Voit & 
Bryan 2001). On the other hand, QSOs/AGNs could supply 
at least part of such a large amount of energy per parti- 
cle, provided that a fraction of their total energy output is 
transferred to the IGM by some physical mechanism, such 
as, for instance, heating by outflows and jets, as suggested by 
several authors (e.g., Valageas & Silk 1999; Inoue & Sasaki 
2001; Yamada & Fujita 2001). Interactions between radio 
lobes of powerful radiosources and the X-ray emitting IGM 
in clusters have been recently studied with Chandra (Fabian 
et al. 2000; Fabian et al. 2001). 



3 THE MODEL 

3.1 Basic assumptions 

Since we are mostly interested in understanding spheroidal 
galaxies with substantial mass {Msph k, Mq) formed 
at high redshift, we concentrate on star formation inside 
laxge virialized haloes of dark matter (DM). At the begin- 
ning, baryonic and DM share the same density profile. Then, 
baryonic matter cools, collapses and starts forming stars. A 
single zone ISM with instantaneous mixing of gas is assumed 
through the inner region where the processes of star for- 
mation from cold gaseous clouds, feedback from SNe, and 
accretion of mass by infall from the surrounding, cooling 
medium are taking place. The star formation is a very effi- 
cient process, building up the bulk of the stellar population 
on a very short time-scale, until the QSO shines at the cen- 
tre. At that time the star formation stops and the galaxy 
undergoes a subsequent phase of passive evolution - Type 
la SNe, exploding in an already ionized medium since then, 
are very efficient in keeping hot the gas (Recchi, Matteucci 
& D'Ercole 2001). The age of the galaxy models, Tgai, is 
a function of the adopted cosmology and redshift of galaxy 
formation, Zf. 

Following Navarro, Frenk & White (1997) and the gen- 
eralization by Bullock et al. (2001), a virialized dark halo of 
mass Mvir identified at z — zo has a virial radius defined by 



47r 3 

Mvir = —Ayir{zo)pu{zo)r^i^ 



(1) 



where pu{zo) is the mean universal density and Avir{zo) is 
the virial overdensity at that redshift. For flat cosmologies a 
good approximation has been derived by Bryan & Norman 
(1998): 



18n^ + 82x - 39x^ 



(2) 



where x = i}{z) — 1 and Q{z) is the ratio of the mean matter 
density to the critical density at redshift z. 

We can also define a virial velocity, V^i^ = GMvir/rvir, 
and the temperature of the gas in hydrostatic equilibrium, 

1 2 

kT = -/imp Kir, (3) 

where p^rrip is the mean molecular weight of the gas. 

Following numerical experiments (see, e.g., Navarro, 
Frenk & White 1996, 1997; Bullock et al. 2001), the den- 
sity profile within a virialized halo is well represented by: 



p{r) 



Ps 



cx{l + cx) 



(4) 



where x = r/rvir and c is the concentration parameter. In- 
tegrating the density profile over the radius up to the virial 
radius we get the halo mass: 



A 3-3 



ln(l + c) 



l+c 



(5) 



Numerical simulations (Navarro et al. 1997; Bullock et al. 
2001) show that the concentration parameter c is a function 
of the mass (c oc M~^'^) and of the redshift [c cx (1 -I- .zo)"^]- 
These dependences have been included in the model. 

In this model we assume that the baryon component in 
virialized haloes has initially the same distribution as the 
DM. However, baryons cool down in a characteristic time 
defined at each radius through the ratio of specific energy 
content to cooling rate, 



, ,^_'ipgas(^_kr[__ 
tcooi(r) - 2 ni(r)A(T)' 



(6) 



where Pgas{r) is the gas density, ne(r) is the electron density, 
nip is the proton mass, piUp is the mean molecular weight 
of the gas and A(T) is a suitable cooling function. In the 
following we adopt cooling functions taken from Sutherland 
& Dopita (1993), which include the dependence on metal 
abundance. 

The second relevant time is the dynamical time, tdyn{r) 
= [3n /32Gp{r)Y'''^ . We neglect the angular momentum of 
the baryons. Actually, the linear tidal torque theory shows 
that later collapsing shells had more time to gain spin via 
tidal torques. As a result, the distribution of the spin param- 
eter exhibits a log normal distribution with the mean value 
of the spin parameter A and its standard deviation ax signif- 
icantly decreasing with increasing redshift (see, e.g.. Mailer, 
Dekel & Somerville 2002). Therefore, a large fraction of low- 
spin objects are expected among the haloes virialized at high 
redshift, z ^ 2-3, which in our model are the hosts of the 
present day massive spheroids (see below) . 

The third time, Tturst = tQSO — tvir, is the time delay 
between the virialization time t^ir, which we assume to be 
also the time t / when stars begin to form in the galaxy, and 
the time of the shining of the QSO tqso, which sets the end 
of the formation of the bulk of the stars. The time Tturst 
is directly related to the visibility of the proto-ellipticals in 
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the far-IR. As shown by Granato et al. (2001), the 850 /^m 
source counts down to 0.5 mjy and related statistics are well 
reproduced by assuming that: 



Tburst(Msph) — 



n, 

t; + log ■ 



HMsph > M, 



sph 
sph 



(7) 



where ~ 0.5 Gyr is the burst time-scale for galaxies with 
Msph > M*p^, and M*^^ ~ 1.5 x lO" Mq is the character- 
istic mass in the local stellar mass function of galaxies. 

For each halo of mass Mvir virialized at redshift z^ir we 
can define a radius rcoii, as the radius within which 



max[tcool{r),tdyn{r)] < Tburst- 



(8) 



The radius rcoii coincides with r^ir for small mass haloes, 
while it is significantly smaller for larger haloes. This is due 
to the increase of the cooling time with increasing halo mass 
and produces a natural cut-off in the stellar mass in galaxies. 
We can also define the time-scale of the gas infall, ri„/ 

= maxltcool{rcoll), tdyn{rcoll)]- 

The fundamental equations of chemical evolution can 
be written as: 

i^Mi).-x,,M,. (t) 

V / inf 



(9) 




where 



R^it) = 



^{t - TM)QAu{t - TM)(t>{M)dM+ 



A 



Mb„ 



4>(Mb) 




(1-A) 



Mb„ 



i^it - TM)QMi{t - TM)<j}{M)dM+ 



ip{t - rM)QAu{t - rM)(t>{M)dM. 



(10) 



Mb„ 



Gi{t) = Xi{t) Mcoid(t) is the cold gas mass in form of the 
element i. The quantity Xi{t) represents the abundance by 
mass of the element i (by definition, the summation over all 
the elements present in the gas mixture is equal to unity). 
The various integrals of Eq.(lO) represent the rates at which 
SNe (I and II) as well as single low- and intermediate-mass 
stars restore their processed and unprocessed material to the 
ISM (see Matteucci & Greggio 1986 for details). 
The star formation rate (SFR) is given by: 



m = 



inayL[tcooi{r),tdyn{r)] 



dr 



dr, 



(11) 



where Mcoid{t) is the gas mass which is cold by time t and 
k = 1. We assume that the gas distribution follows the DM 
distribution during the burst of star formation. Of course, 
the above formula can be rearranged in the form i/)(f) = 



M^oid{t)/T-<. = P-^M'^^i't), with r* being the appropriate 
time-scale for star formation and v-^ the corresponding star 
formation efficiency. 

The infall term is expressed as: 



AGi\ Mgas 1 

at / Tinf 1— e <" 

V / inf 



(12) 



where {Xi)inf — {Xi)p, i.e., the infalling gas has a primor- 
dial chemical composition, and Mgas is the gas mass within 
Tcoii- In this formulation the accretion lasts during all the 
burst, before the QSO shines. 

The last term of Eq.(9) is the rate of reheating, which 
accounts for feedback from SNe. It gives the amount of cold 
gas which is heated and subtracted to further stellar pro- 
cessing per unit time due to SN explosions: 




= X^{t)i:{t)e 



4 rjsNEsN 

5 ~~V^- 



(13) 



(Kauffmann, White & Guiderdoni 1993). risN is the number 
of Type II SNe expected per solar mass of stars formed and 
is computed according to the particular choice of the IMF, 
Es N is the kinetic energy of the ejecta from each supernova 
(~ 10""^ erg) and e is the fraction of this energy which is 
used to reheat the cold gas to the virial temperature of the 
halo. It is worth noticing that we neglect the contribution of 
Type la SNe to the reheating. In fact, during the burst the 
main contribution in terms of energy is coming from Type 
II SNe, and the end of the burst is set by the onset of the 
QSO activity. Thus, we expect that most of Type la SNe 
explode in a warm, rarefied medium. Therefore, they play 
their major role in polluting the intracluster medium (ICM) 
with their end products at later times, rather than efficiently 
regulate the process of star formation at former times. 

An important parameter entering all chemical evolu- 
tion models is the IMF slope, which describes the relative 
numbers of stars born at a given mass. In spite of recent 
observational progress, many fundamental properties of the 
IMF are still unknown. In particular, there is no consensus 
on the question of whether the IMF is independent of en- 
vironmental conditions such as stellar density and metallic- 
ity. New investigations have come to different conclusions on 
the universality of the IMF (Eisenhauer 2001; Gilmore 2001; 
Kroupa 2001); nevertheless it seems to exist some empirical 
indication that the IMF is systematically biased towards 
more massive stars in the early Universe and in starbursts 
(Eisenhauer 2001 and references therein). We will consider 
three cases: a Salpeter (1955) IMF, <j){M) oc M''^-^^ , ii) 
an IMF slightly biased towards more massive stars, (j){M) oc 
M-l■^^ and Hi) a two slope IMF, <^(M) oc M-°-^ for M < 
1 Mq and <j){M) oc M~'^-^^ for M > 1 Mq, aU normalized 
to 1 over the mass range 0.1-100 Mq. Besides, we will as- 
sume that the IMF is not varying with time, since this seems 
to be the case for our own galaxy (Chiappini, Matteucci & 
Padoan 2000). 



3.2 Stellar nucleosynthesis prescriptions 

We compute the abundance evolution of the following ele- 
mental species: H, D, ^He, ''He, '^C, '^C, "N, '"^O, Ne, Mg, 
Si, Fe, and neutron-rich isotopes synthesized from ^^C, ^''C, 
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Table 1. Mg and Fe isotopes in the ejecta of core-collapse supernovae (Nomoto et al. 1997). 



Mi, 



'Mg 



MS, 



13 


9.23 X 10" 


'3 


1.38 


X 


10" 


-3 


8.96 


X 


10" 


-4 


2.10 


X 


10" 


-3 


1.50 


X 


10- 


-1 


4.86 


X 


10- 


-3 


3.93 


X 


10" 


9 


15 


3.16 X 10" 


-2 


2.55 


X 


10" 


-3 


2.03 


X 


10" 


^3 


4.49 


X 


10" 


-3 


1.44 


X 


10- 


-1 


4.90 


X 


10- 


-3 


1.27 


X 


10" 


■8 


18 


3.62 X 10" 


-2 


7.54 


X 


10" 


-3 


5.94 


X 


10" 


-3 


6.04 


X 


10" 


-3 


7.57 


X 


10" 


-2 


2.17 


X 


10" 


-3 


1.37 


X 


10" 


■8 


20 


1.47 X 10" 


-1 


1.85 


X 


10" 


-2 


1.74 


X 


10" 


-2 


2.52 


X 


10" 


-3 


7.32 


X 


10" 


-2 


3.07 


X 


10" 


3 


3.70 


X 


10" 


9 


25 


1.59 X 10" 


-1 


3.92 


X 


10" 


-2 


3.17 


X 


10" 


-2 


4.81 


X 


10" 


'3 


5.24 


X 


10" 


-2 


1.16 


X 


10- 


'3 


8.34 


X 


10- 


9 


40 


3.54 X 10" 


-1 


4.81 


X 


10" 


^2 


1.07 


X 


10" 


-1 


9.17 


X 


10" 


-3 


7.50 


X 


10- 


-2 


2.29 


X 


IQ- 


-3 


1.29 


X 


IQ- 


■8 


70 


7.87 X 10" 


-1 


1.01 


X 


10- 


^1 


2.91 


X 


10" 


-1 


5.81 


X 


10- 


-3 


7.50 


X 


10- 


-2 


3.83 


X 


10- 


-3 


4.17 


X 


10- 


■8 


100 


1.22 




1.50 


X 


10" 


'1 


4.80 


X 


10" 


-1 


5.81 


X 


10" 


-3 


7.50 


X 


10" 


-2 


3.83 


X 


10- 


'3 


4.17 


X 


10" 


■8 



Table 2. Model parameters: dark halo mass, baryonic mass inside the halo, fraction of the baryonic mass which cools 
and collapses, fraction of the kinetic energy from Type II SNe which reheates the gaseous medium, circular velocity of 
the halo at the time of galaxy formation, gas temperature at the time of galaxy formation, efficiency of star formation, 
duration of the burst of star formation. The parameters listed here are relevant to the case zj = 5. The masses are 
expressed in units of Mq, V^ir is in units of km s~^, T is in units of K, Vi, is in units of Gyr~^, and Ti,^,.^^ = tqso ~ 
is in units of Gyr. See text for more details. 



Model 




M 


Mbar 


a 




e 






T 






^burst 


la 


1.65 


X 


1013 


2.47 


X 


1012 


0.21 


0, 


,15 


749.4 


2.04 


X 


lO'' 


5.94 


0.60 


2a 


2.00 


X 


1012 


3.00 


X 


10" 


0.77 


0, 


.15 


370.9 


4.99 


X 


10^ 


7.27 


0.86 


3a 


6.00 


X 


10" 


9.00 


X 


lOlO 


1.00 


0, 


.15 


248.3 


2.24 


X 


10^ 


8.40 


1.17 


4a 


3.00 


X 


10" 


4.50 


X 


lOlO 


1.00 


0, 


.15 


197.1 


1.41 


X 


10^ 


8.63 


1.56 


5a 


2.00 


X 


10" 


3.00 


X 


lOlO 


1.00 


0, 


,15 


172.2 


1.08 


X 


10^ 


8.77 


1.86 


6a 


1.26 


X 


1 nil 

10 


1.89 


X 


10 


1.00 


0, 


,15 


147.6 


7.91 


X 


10® 


8.93 


2.09 


7a 


7.50 


X 


lOlo 


1.13 


X 


lOlO 


1.00 


0, 


,15 


124.2 


5.59 


X 


10^5 


9.13 


2.47 


8a 


2.10 


X 


10"' 


3.15 


X 


109 


1.00 


0, 


,15 


81.2 


2.39 


X 


10® 


9.65 


3.32 


16 


1.59 


X 


1013 


2.39 


X 


1012 


0.21 


0, 


,10 


740.2 


1.99 


X 


10^ 


6.00 


0.60 


26 


1.80 


X 


1012 


2.70 


X 


10" 


0.80 


0, 


,10 


358.1 


4.65 


X 


10^ 


7.32 


0.86 


36 


5.25 


X 


10" 


7.87 


X 


lOlO 


1.00 


0, 


,10 


237.5 


2.05 


X 


lO'' 


8.44 


1.17 


46 


2.60 


X 


10" 


3.90 


X 


lOlO 


1.00 


0, 


,10 


187.9 


1.28 


X 


10*^ 


8.68 


1.56 


56 


1.70 


X 


10" 


2.55 


X 


lOlO 


1.00 


0, 


,10 


163.1 


9.65 


X 


10® 


8.82 


1.86 


66 


1.05 


X 


10" 


1.57 


X 


lOio 


1.00 


0, 


,10 


138.9 


7.00 


X 


10® 


9.00 


2.09 


76 


6.13 


X 


lOio 


9.19 


X 


109 


1.00 


0, 


,10 


116.1 


4.89 


X 


10® 


9.21 


2.47 


86 


1.68 


X 


lOio 


2.52 


X 


109 


1.00 


0, 


,10 


75.4 


2.06 


X 


10® 


9.75 


3.32 


Ic 


1.50 


X 


1013 


2.25 


X 


1012 


0.22 


0, 


,10 


726.0 


1.91 


X 


10^ 


6.07 


0.60 


2c 


2.00 


X 


1012 


3.00 


X 


10" 


0.77 


0, 


,10 


370.9 


4.99 


X 


lO'^ 


7.27 


0.86 


3c 


6.45 


X 


10" 


9.67 


X 


lOlO 


1.00 


0, 


,10 


254.4 


2.35 


X 


10^ 


8.34 


1.17 


4c 


3.28 


X 


10" 


4.92 


X 


lOio 


1.00 


0, 


,10 


203.0 


1.50 


X 


10^ 


8.60 


1.56 


5c 


2.22 


X 


10" 


3.33 


X 


lOio 


1.00 


0, 


,10 


178.3 


1.15 


X 


10*^ 


8.73 


1.86 


6c 


1.42 


X 


10" 


2.12 


X 


lOio 


1.00 


0, 


,10 


153.5 


8.55 


X 


10® 


8.89 


2.09 


7c 


8.50 


X 


lOio 


1.27 


X 


lOio 


1.00 


0, 


,10 


129.5 


6.08 


X 


10® 


9.08 


2.47 


8c 


2.45 


X 


lOio 


3.67 


X 


10^ 


1.00 


0, 


,10 


85.5 


2.65 


X 


10® 


9.59 


3.32 


Id 


1.37 


X 


1013 


2.05 


X 


1012 


0.24 


0, 


,10 


703.5 


1.80 


X 


107 


6.19 


0.60 


2d 


1.90 


X 


1012 


2.85 


X 


10" 


0.78 


0, 


,10 


364.6 


4.83 


X 


10*^ 


7.30 


0.86 


3(i 


6.20 


X 


10" 


9.30 


X 


lOio 


1.00 


0, 


,10 


251.0 


2.29 


X 


10** 


8.36 


1.17 


Ad 


3.16 


X 


10" 


4.74 


X 


lO" 


1.00 


0, 


,10 


200.5 


1.46 


X 


10^ 


8.61 


1.56 


5d 


2.15 


X 


10" 


3.23 


X 


lOio 


1.00 


0, 


,10 


176.4 


1.13 


X 


10*^ 


8.74 


1.86 


6d 


1.37 


X 


10" 


2.05 


X 


lOio 


1.00 


0, 


,10 


151.7 


8.35 


X 




8.90 


2.09 


7d 


8.25 


X 


lOio 


1.24 


X 


lOio 


1.00 


0, 


,10 


128.2 


5.96 


X 


10® 


9.09 


2.47 


8d 


2.39 


X 


lOio 


3.59 


X 


109 


1.00 


0, 


,10 


84.9 


2.61 


X 


10® 


9.60 


3.32 



N, and O. Starting from an initial primordial chemical 
composition of 24 per cent ''He and the remaining hydrogen, 
the evolution of each element is followed by using the for- 
malism of the production matrix Qij{M), first introduced by 
Talbot & Arnett (1973), that gives the fraction of the mass 
of an element j initially present in a star of mass M that 



is transformed in element i and ejected. Detailed nucleosyn- 
thesis prescriptions are from: i) van den Hoek & Groenewe- 
gen (1997) for low- and intermediate- mass stars (0.8-8 Mq) 
(their case with variable mass loss scaling parameter tjagb 
for stars on the asymptotic giant branch); ii) Charbonnel 
& do Nascimento (1998) for 3He production/destruction in 
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Figure 1. Masses ejected in form of ^*Mg (solid squares and 
continuous line) and in form of ^*Mg + -^^Mg + ^®Mg (open 
squares and dashed line) as a function of the initial mass of the 
progenitor star (the values relevant to the 100 Mq star have been 
extrapolated). Yields from Nomoto et al. (1997). 



0,15 



0.1 

CD tC, 



0.05 




Figure 2. Masses ejected in form of ^^Fe (solid squares and con- 
tinuous line) and in form of ^*Fe + ^^Fe + ^''Fe + ^^Fe (open 
squares and dashed line) as a function of the initial mass of the 
progenitor star (the values relevant to the 100 Mq star have been 
extrapolated). Yields from Nomoto et al. (1997). 



low-mass stars (Al < 2 Mq); tii) Nomoto et al. (1997) for 
Type II SNe (M > 13 Mq); iv) Thielemann, Nomoto & 
Hashimoto (1993) for Type la SNe (exploding white dwarfs 
in binary systems). 

There has been some claim in the literature that in order 
to compare theoretical resuhs with observations one should 
consider all isotopes rather than only the dominant ones as 
usually done in chemical evolution studies (see, e.g., Gibson 
& Woolaston 1998) . In Fig. 1 we compare the magnesium 
masses ejected in form of the main isotope with those ejected 
in form of all isotopes, as a function of the initial mass of the 
progenitor star. In Fig. 2 we do the same for iron. The same 
quantities are listed in Table 1. Notice that values for a 100 
Mq star have been extrapolated. As far as iron is concerned, 
it is nearly irrelevant to add the less abundant isotopes to 
the main one, whereas this is not true for magnesium. We 
computed both models including all isotopes and models 
including only the main ones. 



3.3 The average metallicity of a composite stellar 
population 

The average metallicity (or abundance, in general) which 
should be compared with the indices should be averaged on 
the visual hght, namely: 



i /light — 



(14) 



Xi and luminosity Lv ■ On the other hand, the real average 
abundance should be the mass-averaged one, namely: 



{Xi)rnass — ^ 

Otot 



Sto 



X^{S)dS 



(15) 



where Stot is the total mass of stars ever born (Pagel & 
Patchett 1975). Here we will mostly use this last equation, 
since for massive elliptical galaxies the difference between 
the mass-averaged metallicity and the luminosity- averaged 
one is almost negligible (Yoshii & Arimoto 1987; Gibson 
1997). 



4 RESULTS 

4.1 The chemistry of the gas and the stars 

In the reference case SIa ~ 0.7, = 0.3, h — 0.7, Zf — 
5, four sets of eight models each have been computed (it 
is worth noticing that Zf can span a large redshift inter- 
val, the corresponding volume density of virialized haloes 
being specified by the Press- Schechter formalism - Press & 
Schechter 1974; Sheth & Tormen 1999). The model parame- 
ters are shown in Table 2. The first two sets (labelled a and 
b) are relevant to a Salpeter IMF, (t>{M) oc M"^'^'^. Set c as- 
sumes a flatter IMF, (f){M) cx M~^ -^^, and set d refers to a 
two slope IMF, <^(M) oc M'"'^ for M < 1 Mq and (/)(A'/) oc 
j^,j--i.25 £qj. y ^ Mq. Set a refers to the case in which the 
efficiency of reheating from SNell, e, is 15 per cent, whereas 
sets b, c, and d refer to a case with a slightly lower efficiency 
of reheating from SNell, namely 10 per cent. 

Before discussing model results, we want briefiy to fo- 
cus on the behaviour of the star formation efficiency, i^* , as a 
function of the galactic mass, i^* has been computed accord- 
ing to Eq.(ll). Ui, turns out to be a decreasing function of 
the galactic mass (see Table 2, column 8). Nevertheless, our 
star formation efficiencies yield average star formation rates 
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Table 3. Model results (case zj = 5). Mi,{tQso) is the stellar mass that is formed until Iqso is reached, expressed in 
units of Mq. The correction factors 1.6 (Salpeter IMF) and 2.4 (two slope IMF) should be applied in order to obtain 
the mass in stars dead and alive at the present time (Mspft).The chemical properties of the composite stellar population 
are averaged on the mass. The metallicity indices Mg2 and (Fe) given in the last two columns are computed following 
the prescriptions of Matteucci et al. (1998). See text for more details. 



ivioaei 


"'i*(tQ.SO) 


/ rT7„ /Tji \ 
([lo/H]) 


([Mg/H]) 


([Mg/lo]) 


/FT? /XJ1\ 


/FT? /T?^l\ 

([t/lo]) 






Mg2 


{he) 


la 


2.92 X 10" 


-0.420 


-0.038 


0.382 


-0.113 


0.307 


-0.103 


0.0137 


0.232 


2.563 


2a 


1.24 X 10" 


-0.343 


-0.018 


0.325 


-0.083 


0.260 


-0.072 


0.0147 


0.238 


2.656 


3a 


5.94 X 10"' 


-0.345 


-0.084 


0.261 


-0.140 


0.205 


-0.126 


0.0131 


0.226 


2.624 


4a 


2.44 X lO^o 


-0.396 


-0.177 


0.219 


-0.226 


0.170 


-0.210 


0.0108 


0.209 


2.525 


5a 


1.40 X 10^° 


-0.442 


-0.245 


0.197 


-0.290 


0.151 


-0.274 


0.0094 


0.198 


2.441 


6a 


7.21 X 10^ 


-0.517 


-0.333 


0.184 


-0.377 


0.140 


-0.359 


0.0077 


0.182 


2.315 


7a 


3.34 X 10^ 


-0.609 


-0.443 


0.166 


-0.484 


0.125 


-0.466 


0.0061 


0.164 


2.158 


8a 


4.61 X 10** 


-0.888 


-0.750 


0.138 


-0.787 


0.101 


-0.767 


0.0031 


0.112 


1.693 


16 


2.97 X 10" 


-0.413 


-0.031 


0.382 


-0.105 


0.308 


-0.096 


0.0140 


0.234 


2.574 


26 


1.24 X 10" 


-0.317 


0.010 


0.327 


-0.056 


0.261 


-0.045 


0.0156 


0.243 


2.697 


36 


5.90 X 10^" 


-0.293 


-0.029 


0.263 


-0.085 


0.207 


-0.072 


0.0147 


0.236 


2.707 


46 


2.49 X lO^'' 


-0.326 


-0.106 


0.221 


-0.155 


0.171 


-0.140 


0.0127 


0.222 


2.635 


56 


1.42 X 10^° 


-0.365 


-0.167 


0.198 


-0.213 


0.152 


-0.197 


0.0111 


0.211 


2.563 


66 


7.27 X 10^ 


-0.434 


-0.250 


0.184 


-0.294 


0.141 


-0.276 


0.0093 


0.197 


2.447 


76 


3.34 X lO'' 


-0.521 


-0.355 


0.166 


-0.396 


0.125 


-0.378 


0.0074 


0.179 


2.299 


86 


4.62 X 10** 


-0.791 


-0.653 


0.138 


-0.689 


0.101 


-0.669 


0.0038 


0.129 


1.851 


Ic 


3.04 X 10" 


-0.195 


0.288 


0.483 


0.209 


0.405 


0.215 


0.0275 


0.297 


2.949 


2c 


1.26 X 10" 


-0.153 


0.289 


0.442 


0.215 


0.369 


0.222 


0.0279 


0.297 


2.996 


3c 


5.95 X lO^" 


-0.201 


0.197 


0.398 


0.129 


0.330 


0.137 


0.0233 


0.279 


2.906 


4c 


2.48 X lO^'' 


-0.275 


0.089 


0.364 


0.026 


0.301 


0.036 


0.0187 


0.258 


2.779 


5c 


1.42 X 10^° 


-0.335 


0.013 


0.348 


-0.048 


0.287 


-0.037 


0.0159 


0.243 


2.679 


6c 


7.25 X lO'' 


-0.423 


-0.084 


0.339 


-0.144 


0.279 


-0.133 


0.0129 


0.225 


2.539 


7c 


3.32 X 10^ 


-0.528 


-0.202 


0.326 


-0.261 


0.267 


-0.249 


0.0099 


0.202 


2.368 


8c 


4.62 X 10** 


-0.825 


-0.520 


0.305 


-0.576 


0.249 


-0.563 


0.0049 


0.143 


1.895 


Id 


3.04 X 10" 


-0.132 


0.327 


0.459 


0.250 


0.383 


0.256 


0.0299 


0.305 


3.034 


2d 


1.26 X 10" 


-0.093 


0.329 


0.423 


0.260 


0.353 


0.266 


0.0305 


0.306 


3.080 


3d 


5.98 X 10^" 


-0.143 


0.240 


0.383 


0.177 


0.320 


0.185 


0.0257 


0.288 


2.989 


4d 


2.48 X 10^" 


-0.220 


0.132 


0.353 


0.075 


0.295 


0.084 


0.0207 


0.266 


2.85 


5d 


1.42 X 10^" 


-0.282 


0.056 


0.338 


0.001 


0.283 


0.011 


0.0177 


0.252 


2.756 


6d 


7.23 X 10^ 


-0.372 


-0.044 


0.329 


-0.097 


0.275 


-0.086 


0.0143 


0.233 


2.612 


7d 


3.32 X 10^ 


-0.479 


-0.163 


0.317 


-0.214 


0.265 


-0.203 


0.0110 


0.210 


2.439 


8d 


4.64 X 10** 


-0.779 


-0.482 


0.297 


-0.530 


0.248 


-0.518 


0.0054 


0.151 


1.962 



Table 4. Same as Table 3 in the case of Zf = 9, for set d. 



Model 


M,(tQso) 


([Fc/H]) 


([Mg/H]) 


([Mg/Fc]) 


([E/H]) 


([E/Fe]) 


([Z/H]> 


(Z) 


Mg2 


(Fe) 


Id 


3.05 X 10" 


-0.136 


0.345 


0.481 


0.265 


0.401 


0.269 


0.0308 


0.309 


3.037 


2d 


1.26 X 10" 


-0.025 


0.396 


0.421 


0.327 


0.352 


0.333 


0.0351 


0.319 


3.183 


3d 


5.96 X lOif 


-0.032 


0.348 


0.380 


0.287 


0.319 


0.294 


0.0324 


0.309 


3.159 


id 


2.49 X lO^" 


-0.100 


0.255 


0.355 


0.198 


0.298 


0.207 


0.0269 


0.291 


3.045 


bd 


1.43 X 10"' 


-0.156 


0.187 


0.343 


0.132 


0.288 


0.142 


0.0234 


0.277 


2.953 


6d 


7.24 X 10^ 


-0.241 


0.095 


0.335 


0.042 


0.283 


0.052 


0.0193 


0.259 


2.819 


7d 


3.32 X 10^ 


-0.344 


-0.018 


0.325 


-0.069 


0.274 


-0.059 


0.0151 


0.237 


2.655 


8d 


4.64 X 10** 


-0.636 


-0.327 


0.309 


-0.376 


0.261 


-0.364 


0.0077 


0.179 


2.191 



increasing with increasing the galactic mass - (SFR) oc Ml^y^ 
- and guarantee the correct trend of increasing [Mg/Fe] ra- 
tios with increasing total galactic mass (see Table 3). This 
latter result is mainly due to the implementation in Eq.(9) of 
a reheating term proportional to the rate of star formation, 
which prevents smaller galaxies from efficiently converting 
gas into stellar generations at too early times. Owing to that, 



in smaller galaxies a significant fraction of stars form later 
on, out of gas significantly depleted in Type II SN products 
and enriched in iron by Type la SN explosions. Conversely, 
the short duration of the burst secures a lower relative abun- 
dance of iron with respect to magnesium in larger galaxies. 

In Fig. 3 we show the star formation rate, the stellar 
mass and the metal enrichment of both the ISM and the 



© 2002 RAS, MNRAS 000, 0-?? 



Chemical evolution and joint formation of quasars and spheroids 9 




7 8 9 10 543210 

log T (yr) ^ 



Figure 3. Star formation rate, cumulative mass in stars, metal enrichment of gas (lines) and stellar component (stars) 
as functions of time from the beginning of galaxy formation (left panels) and redshift (right panels) for Models 8b, 5b, 
and 16 (M^pjj ~ 3 X 10* Mq, 10^^ Mq, and 2 X 10^^ Mq, respectively, from bottom to top in each panel, except for 
panels e) and f) where only results for Models 8b and lb are shown). The arrows mark the QSO shining time/redshift. 
The adopted IMF is the Salpeter one. 



stellar component as functions of time and redshift for some 
selected models. The models in the figure refer to present 
day galaxies with masses ranging from ~ 3 x 10^ Mq to 
~ 2 X W^^Mq, and to the Salpeter IMF case. The formation 
of early-type galaxies proceeds as a maximum intensity star- 
burst (see also Elmegreen 1999) , reaching a rate as high as ~ 
1000 Mq yr"^ in the case of the most massive spheroids [up- 
per lines in panels a) and b) of Fig. 3], until the QSO shines 
at the centre. At that time, the star formation ceases. The 
chemical properties of the stellar populations, listed in Ta- 
ble 3 for the enlarged grid of masses, have been computed at 
the present time, assuming passive evolution since t = tqso- 
It can be seen that, while by adopting a Salpeter IMF one 
can hardly attain in the gas the super-solar metallicity in- 
ferred from QSO spectra, flattening the IMF is an easy way 
to overcome this problem (see also Matteucci & Padovani 



1993). Results relevant to sets b (Salpeter IMF) and d (two 
slope IMF) are displayed in Fig. 4: at the time of the shin- 
ing of the QSO (indicated by the arrow), Zgas is Zq for 
Model lb {Msph - 2 x 10"Mq) and ~ 3 for Model 
Id {Msph ~ IO^^Mq). These values are in good agreement 
with the range inferred from the analysis of the most ro- 
bust diagnostics in QSO spectra, such as N III] / O III] and 
NV/(ClV+OVl) (Hamann et al. 2001). The actual QSO 
metallicities could be as much as 2 - 3 times higher than the 
above estimates (Hamann et al. 2001; Warner et al. 2001). 
However, it should be stressed that observed values refer to 
the BELRs and reflect the composition of the very central 
region of the host galaxy. Our predictions, derived with a 
one-zone chemical evolution model, are indeed lower limits 
for the central regions. 

As far as the low-luminosity galaxy is concerned (Mod- 
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Figure 4. Metal enrichment of gas (lines) and stars (symbols) as functions of time from the beginning of galaxy 
formation (left panels) and redshift (right panels) for Models 8b (M^ph ~ 3 X 10* A/©) and 16 {Msph ~ 2 X 10^'^ Mq) 
[from bottom to top, panels a) and b)] and Models 8d {M^ph ~ 2 X IO^Mq) and Id {Mapu ~ 1-5 X IO^^Mq) [from 
bottom to top, panels c) and d)]. Set h refers to a Salpeter IMF, </>(Af) oc M~^'^^ , set d to a two slope one of the kind: 
0(M) (X fQj, jv/ < 1 M0 and <I>(M) oc M-^ ^S for M > 1 Mq. The arrows indicate the QSO shining time/redshift 

{tQSO -tf = 0-6 Gyr for Models lb,d; tqso - */ = 3.32 Gyr for Models &h,d). 



els 8 h,d\ Msph ^ 2-3 x 10* M©), negligible differences 
are found by changing the IMF slope. This is due to the 
role of the stellar feedback in low-mass galaxies, which be- 
comes increasingly important when favoring massive stars 
with respect to low- and intermediate-mass stars in the IMF: 
whereas in high-mass spheroids the formation of the stellar 
bulge is ruled by the cooling of the gas and the shining of 
the QSO, in low-mass spheroids the building up of the stellar 
population is rather controlled by the feedback. 

Recent estimates suggest values of [Mg/Fe] between 
0.40 and 0.00 at the centre of ellipticals (Kuntschner 2000; 
Kuntschner et al. 2001; Terlevich & Forbes 2001). These 
values are derived from central absorption-line strengths. 
Corrections for line-strength gradients have to be applied 
in order to get the integrated or composite indices, rep- 



resentative of the mean stellar metallicities, which can be 
compared with results from composite models representing 
whole galaxies. Kobayashi & Arimoto (1999), by using the 
metallicity gradients of 80 elliptical galaxies and by assum- 
ing that Mg2 and Fei reflect the abundances of magnesium 
and iron, respectively, find ([Mg/Fe]) ~ -fO.2 in most of 
early- type galaxies. Their {[Mg/Fe]) does not correlate with 
galaxy mass tracers, at variance with the central [Mg/Fe]. 
However, as the authors themselves stress, the ratio of Mg2 
to Fei may not directly give the [Mg/Fe] ratio. Milone, Bar- 
buy & Schiavon (2000) give instead ([Mg/Fe]) ~ -1-0.3- -1-0.4 
(see their fig. 56). Our theoretical ([Mg/Fe]) are listed in 
Table 3. They display a decreasing trend with decreasing 
galactic mass, even if milder than observed at the galactic 
centre. The ([Mg/Fe]) ratios listed in Table 3 have been 
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Figure 5. [C, N, Mg, Si/Fe] vs. [Fe/H] in tlie gas from t = tf to 
t = tQso foi" Models 16 (continuous lines) and Id (dashed lines). 
Different elements show different degrees of enhancement, owing 
to their different nucleosynthesis histories. 



computed by using the yields of ^"Mg^- Mg+^^Mg and 
''''Fe+'^'^Fe+''''Fe+^**Fe, rather than those of the main iso- 
topes alone, following arguments by Gibson & Woolaston 
(1998). The elemental {[^''Mg/^^Fe]) ratios turn out to be 
only less than 0.10 dex lower on the average. Owing to the 
observational uncertainties, we can not discriminate between 
the two choices. 

If we consider the ratio ([E/Fe]) rather than ([Mg/Fe]), 
where E refers to all 'enhanced' elements, namely C, N, O, 
Ne, Mg, Si, S following Trager et al. (2000a, their model 4 
- see below), we find milder enhancements. This is due to 
the fact that now we are considering also elemental species 
partially (Si, S) if not even almost entirely (C, N) produced 
by low- and intermediate-mass stars restoring their nucle- 
osynthesis products on longer time-scales than SNell (re- 
sponsible for the whole Mg production) . In Fig. 5 we show 
the temporal evolution of the abundance ratios of some en- 
hanced species with respect to iron (from t — tf to t = 
tQSo) for Models 1 b,d. Different enhanced species behave 
very differently; by the way, C should not be included in the 
enhanced group at all. 

Kuntschner & Davies (1998), from a comparison of their 
measurements of central C4668 and H^^ indices in a sam- 
ple of early-type galaxies in the Fornax cluster to predic- 
tions from single-burst stellar population models (Worthey 
1994; Worthey & Ottaviani 1997), derive that their ellipti- 
cals have [Fe/H] from -0.1 to -1-0.6. Kuntschner (2000) finds 
that early-type galaxies brighter than Mb ~ —17 in the For- 
nax cluster form a sequence in metallicity varying roughly 
from —0.25 to -1-0.30 in central [Fe/H]. Making an attempt 
at estimating mean values, Kobayashi & Arimoto (1999) find 
that typical mean stellar metallicities are ([Fe/H]) ~ —0.3 
and range from ([Fe/H]) ~ —0.8 to -1-0.3, well below the 
highest values observed at the galactic centres. Our predic- 



tions (Table 3, third column) are indeed in good agreement 
with the values inferred from the metallicity gradients, es- 
pecially if an IMF fiatter than the Salpeter one is assumed. 
However, we never recover the highest metallicity values. 
Undoubtedly, the fact that in the framework of our model we 
can not assembly objects more massive than Mgph ~ 1.5- 
2 X 10^^ Mq is playing a crucial role (however, it should 
be noticed that the maximum stellar mass which can be 
assembled depends also on the assumed redshift of galaxy 
formation - see Sect. 4. 3 below). 

Very recently, Trager et al. (2000a) have derived single- 
burst stellar population equivalent ages, metallicities, and 
abundance ratios for a sample of local early-type galaxies 
from H/3, Mg&, and (Fe) line strenghts using an extension of 
the Worthey (1994) models that accounts for the enhance- 
ments of Mg and other a-elements relative to the Fe-peak 
elements. The metallicities and enhancement ratios [E/Fe] 
they have found are strongly peaked around -1-0.26 and 
-1-0.20, respectively, in an aperture of radius re/8. Gradients 
in stellar populations within galaxies are found to be mild, 
with metallicity decreasing by 0.20 dex and [E/Fe] remaining 
nearly constant out to an aperture of radius re/2 for nearly 
all systems. In Fig. 6 we compare our ([E/H]) vs. ([Fe/H]) re- 
lation for Models from 16, d to 8b, d (spheroid masses range 
from Msph ~ 2 X IO^Mq to M^ph ~ 2 x 10** M©; set b: 
Salpeter IMF, set d: two slope IMF) to that by Trager et 
al. (2000a, b). Results relevant to Models from Ic to 8c (sin- 
gle slope IMF flatter than the Salpeter) are very similar to 
those relevant to Models from Id to 8d and therefore have 
not been plotted in the graph. As expected, our metallicities 
and abundance ratios agree much better with the Trager et 
al.'s ones computed through the global re/2 aperture rather 
than with those computed through the central re/8 aper- 
ture. Again, the discrepancy is mainly due to the fact that 
our abundances are accounting for the distribution of stel- 
lar populations over the whole galactic physical dimension. 
Besides this, a remarkable point is that the slope of our 
([E/H]) vs. ([Fe/H]) relation is quite similar to that found 
by Trager et al. (2000a,b) for central values. The situation 
is also improving if galaxy formation is pushed at higher 
redshift (Table 4). 

So far, we made an attempt to compare theoretical 
abundances for integrated early-type galaxies with estimates 
coming from the analysis of observed spectral indices. Now 
viceversa we want to use some calibrations and translate 
the theoretical abundance ratios into metallicity indices to 
be directly compared with observations. When dealing with 
early- type galaxies, the most widely used tool to compare 
theoretical predictions on a-element and Fe abundances in 
stars to observations is the Mg2 vs. (Fe) diagram. In Fig. 7 
observations (Worthey et al. 1992; Gonzalez 1993; CaroUo 
& Danziger 1994a, b; Kuntschner et al. 2001) are shown as 
long as model predictions. The theoretical indices have been 
obtained by converting the ([Mg/Fe]) and ([Fe/H]) ratios 
predicted by Models from 16, d to 8b, d listed in Table 3 into 
the metallicity indices Mg2 and (Fe) by means of the cali- 
brations by Tantalo et al. (1998), as described in Matteucci, 
Ponzone & Gibson (1998). The spread in the Mg2 data is 
reproduced remarkably well. In particular, the adoption of a 
two slope IMF of the kind of that described in Sect. 3.1 allows 
us to cover a wider range of values in both indices. This is an 
unprecedent result, particularly interesting owing to the fact 
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Figure 6. The ([E/H]> vs. ([Fc/H]) relations obtained for Models 
from lb,d to 8b,d (M^ph ~ 2 X lO" -2 X IO^Mq; set b: Salpeter 
IMF, set d: two slope IMF) are compared to that found by Trager 
et al. (2000a, b) from H/3, Mgb, and (Fe) measurements in a sample 
of local ellipticals. Our models show a lower Fe and E content, 
owing to the fact that we are running one zone models. Results 
relevant to Models from Ic to 8 c (single slope IMF flatter than 
the Salpeter one) are very similar to those relevant to Models 
from Id to 8d and are therefore not shown here. 
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Figure 7. Mg2 vs. (Fe) theoretical relations compared to the 
available data. The recent, homogeneous sample by Kuntschner 
et al. (2001) is shown as filled circles; older data are shown as open 
circles. Note that the Kuntschner et al. sample refers mostly to 
group and cluster environments. Results for sets b (stars) and d 
(triangles) are shown. The two sets differ uniquely in the IMF 
slope. Results relevant to set c, characterized by a single slope 
IMF flatter than the Salpeter one, are very similar to those rele- 
vant to set d, so we do not show them in the graph. 



that an IMF of the kind adopted here is supported both ob- 
servationally in the low-mass range (e.g., Wyse et al. 1999; 
Zoccali et al. 2000) and theoretically (Elmegreen 2000a,b). 
Note that an IMF with a single slope of 1.15 over the entire 
mass range would produce similar results. The high-mass 
end of the diagram, which is actually missed by our theoret- 
ical models, could be easily covered by adopting an Arimoto 
& Yoshii IMF, cj){M) oc M~°-^^ , i.e., an IMF with an even 
flatter slope. However, as we already noted before, an off- 
set between model predictions and data has to be expected, 
since data refer to central index values, whereas model pre- 
dictions are relevant to the whole physical dimension of the 
galaxy. Milone et al. (2000) find that the integrated indexes 
of their sample galaxies span a range from ~ 0.2 to ~ 0.3 
in Mg2, and from 2 to ~ 2.5 in (Fe), while central index 
values range from 0.25 to ~ 0.35 for Mg2, and from 
2.5 to ~ 3.5 for (Fe). The offset we find is nearly the same 
(cfr. model predictions with data in Fig. 7). Moreover, as we 
already pointed out, in a monolithic picture for galaxy for- 
mation a natural cut-off in the stellar mass of galaxies should 
be expected, owing to the increase of the cooling time with 
increasing the host halo mass. 

Finally, one could inquire about the properties of the 
very rare ellipticals formed at the highest redshifts, which 
are probably the hosts of the highest redshift QSOs {z ^ 6). 
In Table 4 we summarize the results relevant to set d, for 
the case Zf = 9 and = 0.4 Gyr. 



4.2 Broad-band photometry 

We have checked that also the photometric properties of our 
elliptical galaxy models are consistent with observations of 
local galaxies. To this aim, we have used the spectrophoto- 
metric code by Silva et al. (1998) and computed the spectral 
energy distribution (SED) and colours expected according 
to the star formation and chemical evolution histories of the 
models. 

For the shape of the SED, an example is shown in Fig. 8, 
where we compare with the template SED for ellipticals by 
Arimoto (1996). 

As discussed in Section 2.1, elliptical galaxies are ob- 
served to follow a CMR, with redder colours at increasing 
luminosity. The slope of this relation is interpreted as driven 
mostly by metallicity, rather than age, and its tightness as 
a small age dispersion among the bulk of the stellar popula- 
tions. This interpretation rests on the reproduction, by sev- 
eral authors, of a combined set of observational constraints, 
including broad-band magnitudes and colours, spectral in- 
dexes, such as H/3, Mg2, (Fe), etc., and abundance ratios 
(e.g.. Bower et al. 1992; Bressan et al. 1994, 1996; Tantalo 
et al. 1996; Maraston & Thomas 2000). On the other hand, 
there are evidences, from this same data set for local galax- 
ies and from higher- z observations, that at least a minor 
fraction of the stellar populations of elliptical galaxies may 
have formed someway after the bulk of their stars was al- 
ready formed (e.g., Bressan et al. 1996; Franceschini et al. 
1998; Longhetti et al. 2000; Trager et al. 2000a,b). 

In Fig. 9, the colour - magnitude relation V ~ K vs. A4v 
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Figure 8. The SED of Model Id {M^p^ ~ 1.5 X 10"Mq, two 
slope IMF) at 13 Gyr is compared with the template optical-NIR 
SED for cUiptical galaxies by Arimoto (1996). 



and the J — K vs. V — K colours are shown for models of case 
d (two slope IMF, Table 3) with mass M^ph > 10^°Mq at 
the age of 10 and 15 Gyr (zf = 5 implies an age of 12 Gyr for 
the assumed cosmological parameters). Lower mass models, 
with prolonged star formations, have more prolonged evo- 
lutionary histories, for which we do not have enough con- 
straints in the framework of our model. The V — K colour is 
a strong function of the global metallicity of stars. In agree- 
ment with the results we found for the chemical abundances, 
the colours we obtain with the Salpeter IMF are too blue 
because of the very low metallicities reached (cases a and 
b in Table 3). Instead, in order to reproduce the observed 
V—K colours, an average metallicity greater than solar must 
be reached, for the most massive models, which is the case 
for Models d shown in Fig. 9. As remarked in the previous 
section, since the most massive objects in our model have 
M,ph ~ 1.5xlO"M0, we do not cover the highest luminos- 
ity end of the plot. In Fig. 9 we also show the effect, for the 
two less massive objects in the plot at 15 Gyr, of a small 
burst episode involving 5 per cent of the mass, and taking 
place 4 Gyr ago. It is clear that the slope of the predicted 
CMR may fluctuate around the average value determined 
by metallicity because of these uncertainties. 

The spectrophotometric code provides also the blue lu- 
minosities we need in order to compute the rates of Type la 
SNe at the present time. We find an average present SNIa 

0.3 SNu for 2/ = 5 to 



rate varying from {Rate*'sNia) 

{Rate^SNia) ~ 0.2 SNu for Zf = 9, to be compared with the 
observed value of 0.18 ± 0.06 SNu (Cappellaro et al. 1999). 
Although uncertainties in both observations (e.g., empirical 
corrections for the bias in the inner galactic regions) and 
theory (e.g., fraction of stars of a given mass that enters the 
formation of binary systems which will give rise to Type la 
SN explosions; metallicity effects - which we do not account 



Figure 9. Upper line: The colour — magnitude relation V — K 
vs. Mv for Models d (two slope IMF) is shown at the age of 10 
and 15 Gyr. Data are from Bower et al. (1992) (triangles and 
diamonds for Coma and Virgo cluster ellipticals, respectively) 
and Mobasher et al. (1999) (squares. Coma cluster ellipticals). 
The crosses connected by the dotted line are an example of how 
a small star formation episode may affect the V magnitude (see 
text). Lower panel: J — K vs. V — K, same models as above and 
data from Bower et al. (1992). 



for here) prevent us from drawing any firm conclusion, we 
notice that theoretical results agree quite well with observa- 
tions, at least at the 2 a level. 



4.3 The M/Msph Ratio 

While the luminous component of a galaxy is rather well de- 
fined, the associated halo is not easily identifiable as a well 
defined structure, separated from other close haloes. Never- 
theless, luminous galaxies are good tracers of the mass distri- 
bution. Recent statistical studies of the weak lensing associ- 
ated to nearby 'lens' galaxies allow for a first estimate of the 
galaxy-mass correlation function (McKay et al. 2001). These 
authors point out that the galaxy-mass correlation function 
is independent of the environment, confirming that the mass 
within ~ 300 kpc from a central galaxy is clearly asso- 
ciated with the central object. They also find that elliptical 
galaxies in the luminosity range 1.6 < L^i < 17.7 x 10^" 
Lq exhibit an average {{M{R < 260 kpc)/L^/) = 123 ± 14 h 
Mq/Lq. This result holds for dark matter haloes with 
isothermal density profile, while in the case of NFW profile 
the average value is about a half. Assuming a Salpeter IMF, 
our spectrophotometric model yields M/Lz ~ 2.4 Mq/Lq 
for elliptical galaxies with an age Tgai = 13 Gyr while, for the 
two slope IMF we adopted, M/L^ ~ 1.8 Mq/Lq. In the case 
of a NFW profile and Ho = 70 km s~^ Mpc'\ the resuhs 
obtained by McKay et al. (2001) suggest Mhaio/Msph ^ 17 
for a Salpeter IMF and Mhaio/Msph — 23 for the two slope 
IMF adopted in this paper, for galaxies with mass in the 
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Figure 10. M/Mgph ratios vs. Msph- M is tlie dark halo mass, 
Mgpfi is tlie mass in stars and stellar remnants at the present 
time. The ratios depend on both IMF slope [upper panel; sets a, 
b: Salpeter IMF, set c: single slope IMF flatter than the Salpeter, 
set d: two slope IMF; the efficiency of reheating from SNell is set 
to 10 per cent for all cases, but set a where it is 15 per cent) and 
redshift of galaxy formation (lower panel). 



range 6 x 10^° -6 x lO" Mq. From Fig. 10 it is apparent 
that in the mass interval explored the agreement between 
the results of McKay et al. (2001) and the model predic- 
tions is very good. 

The increase of the ratio Mhaio/Msph with decreasing 
mass in the low-mass range is due to the effect of the stellar 
feedback acting on the gas in shallow potential wells. This 
portion of the curve is sensitive to the fraction of the SN 
energy transferred to the gas, and to the IMF. As a result, 
the fraction of the baryons which form stars is strongly de- 
creasing below a few 10^ Mq, in keeping with observations 
showing that in low luminosity /mass galaxies the relative 
fraction of dark matter to stars progressively increases with 
decreasing luminosity (e.g., Persic, Salucci & Stel 1996). 

At large masses the paroxysmal increase of the ratio 
with stellar mass is mainly due to the increase of the cool- 
ing time of the gas in the outer regions of large haloes. The 
fraction of baryons available for cooling and collapsing into 
stars is fixed by imposing Tt,urst < tqso — tvir [cfr. Eq.(8)]. 
The amount of stars formed is limited by the relatively short 
time lag between the virialization and the QSO appearance 
for massive objects. The Mhaio/Msph ratio is not only sen- 
sitive to Tturst, but also to the formation redshift, larger 
masses being more easily assembled at high redshift, since 
at fixed mass the densities of DM and baryons increase with 
Zvir (see Section 2.1). In the case of z-uir = 9, the largest 
mass in stars in our set d with still acceptable M^aio/Msph 
is Maph ~ 4 X 10^^ Mq (Fig. 10). In this framework, more 
massive galaxies are assembled by coalescence, as suggested 
by their often observed distorted morphologies (e.g., double 
cores). On the other hand, the number density of massive 



haloes strongly decreases with increasing redshift. The result 
is the observed exponential decline of the luminosity /mass 
function at large luminosities/masses. The statistical aspects 
implied by these results will be discussed in a forthcoming 
paper. 

It is worth noticing that the values of e, namely the effi- 
ciency of reheating from SNell, are significantly constrained 
by the results of McKay et al. (2001). If e=l is adopted, as 
it seems to be required in order to get, under very favorable 
hypotheses, ~ 1 KeV of energy per particle injected into 
the IGM by SNe (Valageas & Silk 1999; Kravtsov & Yepes 
2000), a minimum value of M/L^ ~ 250 for L^ = 2x 10^° L© 
and for the case d and Z[ = 5 is recovered. This value ex- 
ceeds by a factor of 4 the result of McKay et al. (2001) for 
a NFW density profile. 



5 DISCUSSION AND FINAL REMARKS 

In this work we concentrated mostly on the chemo- 
photometric properties of the stellar populations of early- 
type galaxies. In particular, we tested specific histories for 
the formation of the spheroids, which turned out to provide 
a good fit to crucial observed relations: the mass - met allicity 
relation, the Mg2 vs. (Fe) diagram, the colour - magnitude 
relation V — K ys. My, the J — K vs. V — K colours and 
the Mhaio/Msph ratio. We found that all these properties 
are well reproduced by using the same dependence of the 
burst duration on galactic mass which is requested in order 
to fit the submillimetre source counts (Granato et al. 2001). 
It is worth noticing that the adopted scenario does hold for 
galaxies which are large enough to host a QSO, i.e., Msph 
^ a few 10^ Mq, less massive objects possibly suffering for 
even more complex evolutionary paths than described above. 
Secondary episodes of star formation in massive galaxies at 
z < 1 would produce only a tiny amount of stars; thus, the 
observed evolution of the FP and the observed [O II] emis- 
sion could be explained (Treu et al. 2001) without any need 
of altering the global chemical properties of the stellar pop- 
ulation. 

At variance with most of the previous models of chem- 
ical evolution of spheroidal galaxies, we accounted for the 
stellar feedback through Eq.(13), i.e., we allowed a fraction 
of gas to be heated and subtracted to the ongoing star for- 
mation process at each time. The reheating term acts in 
the direction of reducing the efficiency of star formation in 
less massive spheroids. This reconciles our star formation 
efficiencies, decreasing with increasing galactic mass, with 
the increase of the star formation efficiency with increas- 
ing galactic mass required by other authors and it yields 
(SFR) oc M}pl. The overall result is similar, none the less 
the physical approach is very different. In our model, the 
star formation is governed by the rate at which baryonic gas 
cools and falls into DM haloes and is inhibited by heating 
from SNe and definitively ended by the QSO peak activity 
combined to the SN feedback. 

We find that the ratios {[Fe/H]), {[Mg/H]), {[Mg/Fe]), 
([E/H]), {[E/Fe]), and ([Z/H]) all increase with galactic mass 
(see Table 3) . This is due to the combined effects of the fol- 
lowing aspects of the model: the duration of star formation 
is shorter with increasing galactic mass; ii) the number of 
stellar generations increases with increasing galactic mass 
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Figure 11. Theoretical Type la SN rate for Model Id {M^ph ~ 
1.5 X 10" M©, two slope IMF) as a function of time since the 
beginning of galaxy formation. The bulk of SNIa explosions is 
expected to occur soon after the shining of the QSO, when the 
galaxy shows up as an extremely red object (ERO). 



((SFR) oc Mgp^); Hi) the stellar feedback is more eflective 
in low-mass galaxies (the amount of SNII-enriched material 
retained by low-mass galaxies is lower, due to the shallower 
potential wells). 

If we relaxe the first hypothesis and set Tturst = 2 Gyr, 
constant with mass, the abundance ratios still increase with 
increasing galactic mass but {[Mg/Fe]) and ([E/Fe]), which 
flatten out (see also Kawata 2001). Therefore, a star- forming 
phase which lasts longer in less massive spheroids is strongly 
needed in order to reproduce the observed trend of increas- 
ing [Mg/Fe] in the nuclei of local ellipticals with increasing 
galactic mass. On the contrary, if we suppress the stellar 
feedback (e = 0), the ratios of all the elements with respect 
to hydrogen, ([el/H]), decrease with increasing mass, while 
the increase in {[el/Fe]) is preserved. 

Flattening the IMF with increasing galactic mass is a 
further way of increasing the quantities ([Fe/H]), ([Mg/H]), 
([Mg/Fe]), ([E/H]), ([E/Fe]), and ([Z/H]) with increasing 
galactic mass; however, in the framework of our model, we 
do not need such a flattening in order to increase the overall 
mean metallicity of the stellar population. A careful analysis 
of Tables 3 and 4 actually reveals that the mean metallicity 
of the stellar population starts slightly decreasing again at 
Msph ^ 10^^ Mq (cfr. the el/H abundance ratios for Models 
1 and 2 in both tables). This is due to the increase of the 
cooling time of the gas in the outer regions of the largest 
virialized haloes. This increase limits the amount of stars 
which are born and hence the metal enrichment of the gas. 

In the framework of the adopted model, massive QSO 
hosts {Msph > 2 X 10^° Mq) experience a relatively short 
{Tburst ~ 0.5-1 Gyr), intense (SFR ~ 50-1000 M© yr"^) 
star- forming phase at high redshift [z ^ 3), during which 
they show up as ultraluminous far-IR galaxies [SCUBA 



galaxies). After the shining of the QSO, they appear as 
EROs, for which we predict solar or higher than solar stel- 
lar metallicities. Moreover, we predict that very high rates 
of SNela should be observed in EROs, since in these objects 
the bulk of SNela explode soon after the shining of the QSO 
(see Fig. 11), in a dust-free medium. Therefore, non dusty 
EROs are preferred targets for searching for high-z SNela. 

Intermediate- and low-mass QSO hosts (Msph < 2 x 
10^° M©) are characterized by a more protracted starburst 
phase [Tturst ~ 1-2.5 Gyr) and show up as Lyman break 
galaxies (LBGs) at z J> 3 (see also Granato et al. 2001). The 
curves relevant to Model 56 in the upper panels of Fig. 3 are 
representative of the history of star formation of a typical 
LEG; more generally, our models predict SFRs between ~ 

5 and ~ 50 M© yr"'^ for this kind of objects. These values 
agree quite well with recent estimates of the star forma- 
tion activity in LBGs by Pettini et al. (2001). Moreover, our 
models predict a value of log(0/H) -I- 12 ~ 8.2-9.2 in the 
gas of LBGs at 2: ~ 3, strongly depending on the adopted 
redshift of galaxy formation and hence on the age of the 
stellar population. The lowest log(0/H) -I- 12 value we give 
is associated to a low-mass spheroid (Model 7) which started 
forming stars at ~ 3 and to a stellar population of age 
0.5 x 10* yr; the highest one is associated to a more massive 
spheroid (Model 4) which started forming stars at 2/ 5 
and to a stellar population of age 10® yr. The lowest values 
we find agree well with the estimates by Pettini et al. (2001). 
This is quite reassuring, since their estimates refer to objects 
with virial masses of about 10^" Mq, and Model 7 indeed 
refers to an object with virial mass Mtar ^ 1.2 X lO^^M© 
(see Table 2). A more detailed study of LBGs will be the 
subject of a forthcoming paper. 

The adopted model implies that the ratio Mbh /Msph 
between the final mass of the black hole responsible for the 
QSO activity and the mass of the host galaxy is not evolv- 
ing with redshift, at variance with respect to semi-analytical 
models (Kauffmann & Haehnelt 2000). The most recent 
observations show a modest, if any, variation of this ratio 
with redshift (Kukula et al. 2001; see also the discussion in 
Granato et al. 2001). 

The ratio Mhaio/Msph predicted by the model depends 
significantly on the adopted value of e, the efficiency of 
the transfer of the SN energy output to the ISM. A very 
good agreement with the findings of McKay et al. (2001), 
A'lhaio/Mgph — 20 for galaxies in the mass range 6 x 10^*^ - 

6 X 10^^ Mq, is found under the assumption e ~ 0.1. The 
combination of potential well and feedback makes the ra- 
tio increasing with decreasing mass. The assumption e=l is 
strongly ruled out in our model by the fact that the mini- 
mum ratio Mhaio/Msph would be 4 times larger than that 
found by McKay et al. (2001). 

This last result is also important in the context of the 
preheating of the IGM in galaxy groups and clusters. Ac- 
tually, it has been suggested that with e = 1 and under 
very favorable conditions, SNe can yield as much as ~ 0.5 - 1 
KeV per particle, thus bringing theoretical models of cluster 
formation into agreement with observations. However, this 
would contrast with the Mhaio/Msph ratio found in galaxies. 
The supernovae are likely to be supplemented by some other 
heating mechanism (Balogh, Babul & Patton 1999; Valageas 
& Silk 1999; Wu et al. 2000; Kravtsov & Yepes 2000). In par- 
ticular, Valageas & Silk propose high-redshift preheating of 
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the IGM by radiation from quasars. It is quite natural that 
the mechanism which accounts for the required heating of 
the entire ICM affects also the evolution of its host. 

In the context of the relationship between the QSO and 
the host galaxy, it is also worth noticing that the model 
predicts a metal abundance of the gas at the epoch of the 
QSO appearance in close agreement with the values inferred 
from the BELRs of high-z QSOs. 

In conclusion, the basic assumption of the model, 
namely that massive spheroidal galaxies evolve rapidly af- 
ter the virialization of their haloes till their QSOs shine, is 
reinforced by the match of the model outcomes with the ob- 
served properties. In particular, in this paper we have shown 
that QSOs can play a major role in determing the observed 
chemo-spectrophotometric properties of spheroidal galaxies 
in a wide range of mass. 
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